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ABSTRACT 

The  topic  of  this  work  is  suggested  by  the  Atomic  Energy  Commission'  s 
controlled  fusion  research,    which  involves  an  attempt  to  confine  a  perfectly 
conducting  plasma  in  a  strong  magnetic  field.     Of  several  confinement  schemes 
proposed,   the  cusped  geometry  model  offers  promise  of  stability  in  exchange 
for  some  leakage  through  the  cusps.     It  is  the  purpose  of  this  dissertation  to 
provide  asymptotic  fornaulas  for  the  cusps  of  the  simplest  cusped  geometry, 
which  is  a  spindle-shaped  figure  consisting  of  two  tapered  circular  cones  set 
back  to  back.     This  figure  thus  has  one  circular  knife -edged  cusp  where  the 
cones  are  joined  and  two  identical  axisymmetric  point  cusps.     The  spindle  is 
supposed  to  contain  a  plasma  which  is  considered  a  fluid  at  rest  and  therefore 
at  constant  pressure.     An  external  vacuum  magnetic  field  due  to  a  quadrupole  at 
infinity  balances  the  fluid  pressure.     The  problem  of  finding  the  shape  of  the 
interface  between  the  fluid  and  the  vacuum  is  stated  as  a  differential  equation 
with  two  boundary  conditions.     Garabedian,    by  analytic  continuation  into  the 
four -dimensional  space  of  two  complex  variables,    restates  this  problem  as 
an  integral  equation  in  terms  of  an  analytic  function  g(z).     In  the  real  plane 
the  trace  of  the  interface  is  given  by  z  =  g(z). 

At  the  edge  cusp,    z  =  z  ,    conformal  mapping  is  used  to  obtain  the  con- 
vergent series  expansion  of  g(z)  „ 

2^     °°  2^ 

z  =  g(z)  =  -  z  +  iC(z-z  )    +  2    a   (z-z  ) 

n  =  4 
At  the  point  cusp  z  =  z_  we  have  a  branch  point  which  greatly  complicates 
matters.     We  extend  work  of  Lewy  to  cover  this  singular  case,    obtaining  there- 
by an -asymptotic  expansion 

^  2Sln{-lln{z-z^))     (-|  +  2jfn^)  /(£n(-£n(z-z^)))^' "^ 


+ 


/(£n(-jen(z-ZQ)))  \ 
\{£n(z-z^))^       / 


z  =  g(z)^z+a  i(z-z    )  {.    ,  )+  2  2 

0  0     ]in(z    Zq)        (jen(z-ZQ))^  {lln{z-z^)f  - ._..    .^,  ,  ,  ^ 

valid  in  the  neighborhood  of  a  point  cusp. 
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INTRODUCTION 

The  topic  of  this  work  is  suggested  by  the  Atomic  Energy  Commission'  s 
controlled  fusion  research  [4],    which  involves  an  attempt  to  confine  a  perfectly 
conducting  plasma  in  a  strong  magnetic  field.     Several  confinement  schemes 
have  been  proposed.     Of  these,   the  cusped  geometry  and  picket  fence  models 
[2,  4,  12]  offer  promise  of  a  stable  figure  of  equilibrium  in  exchange  for  a 
certain  amount  of  leakage  through  the  cusps.     Among  the  matters  of  interest 
in  connection  with  these  cusped  naodels  is  the  question  of  the  asymptotic  shape 
of  the  cusps  for  which,    unfortunately,   no  analytic  description  is  available.     It 
is  the  purpose  of  this  dissertation  to  remedy  that  deficiency  by  providing 
asymptotic  formulas  for  the  cusps  of  the  simplest  cusped  geometry,    which  is 
the  spindle-shaped  configuration  of  figure  1. 

This  axisymmetric  configuration  consists  of  two  identical  tapered  cones 
whose  broad  ends  are  joined  to  form  a  circular  knife-edged  cusp,    generally 
called  the  ring,    or  edge,    cusp.     The  two  mirror  image  pointed  ends  of  the 
resultant  spindle-shaped  figure  are  axisymmetric  point  cusps.     This  figure 
is  supposed  to  contain  a  fully  ionized  plasma  which  we  regard  as  a  perfectly 
conducting  fluid  continuum  at  rest,  and  therefore  at  constant  pressure.     Our 
spindle-shaped  figure  results  from  a  balance  between  this  fluid  pressure  and 
the  force  of  an  external  vacuum  magnetic  field  due  to  a  quadrupole  at  infinity. 

The  problem  of  finding  the  shape  of  the  interface  between  the  fluid  and 
the  surrounding  vacuum  is  stated  mathematically  as  a  differential  equation 
with  two  boundary  conditions,    expressed  in  terms  of  a  potential  ip,    analogous 
to  the  stream  function  of  fluid  dynamics.     The  second  boundary  condition  would 
be  superfluous  if  the  boundary  were  given,    but  because  the  boundary  is  unknown, 
we  are  faced  with  a  free  boundary  problem  which  requires  an  extra  boundary 
condition. 


perfectly 

conducting 

plasma 


/  .  ^  magnetic  field. 


Figure  I.         Plasma  Spindle  in  Magnetic  Field. 
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This  problem  is  so  difficult  that  in  spite  of  long  consideration  it 
remains  unsolved.     Not  only  is  the  problem  difficult,    but  its  formulation 
in  terms  of  a  differential  equation  in  0  is  particularly  awkward  because 
the  equation  contains  a  singular  coefficient  1/y  that  becomes  infinite  at  the 
point  cusps. 

We  wish  to  find  asymptotic  expressions  for  the  cusps  because  of 
the  problem  of  leakage  of  the  plasma  through  these  cusps.     To  do  so  we 
shall  utilize  the  work  of  Garabedian  [10]  who,   by  analytic  continuation 
into  the  four -dimensional  space  of  two  complex  variables  has  replaced  the 
differential  equation  for  ii  with  a  more  amenable  integral  equation  for  an 
unknown  analytic  function  g(z).     The  integral  equation  is  obtained  with  the 
aid  of  a  Riemann  function.     In  the  real  plane  the  integral  equation  defines 
i//,    and  the  trace  of  the  spindle  boundary  yields  a  pattern  consisting  of  four 
mirror  image  curves  T  whose  equations  are  given  in  the  form  z  =  g(z). 
In  order  to  satisfy  the  integral  equation  it  is  necessary  to  choose  the  co- 
efficients and  exponents  of  g(z)  in  a  particular  way. 

At  the  edge  cusp  z  =  z     we  invoke  Schwarz'  s  reflection  principle 
to  obtain  a  convergent  series  expansion  of  g(z).     Thus 

3  n 

2       °°  2 

z  =  g(z)  =  -  z  +  i  C(z-z  )     +  S     a    (z-zj 

1  .    n  1 

n  =  4 


which  becomes 


_C  I       3C'  2  ^ 

X  =  - — -  (y-yJ    -^-^(y-yJ    +o((y-y)    ) 

in  terms  of  x  and  y. 

At  the  point  cusp  z  =  z„  the  contour  integral  which  represents  (//  in 
terms  of  the  Riemann  function  is  greatly  complicated  by  the  presence  of  a 
branch  point.     Lewy  [15]  has  treated  some  similar  problems  from  a  different 

-  7  - 


standpoint.     He  gives  a  general  expansion  which  suggests  the  form  of  g(z). 
His  expansion  is  valid  for  geometries  quite  similar  to  ours  but  is  not  quite 
adequate  here  because  of  the  unusually  singular  nature  of  our  point  cusp. 
Nevertheless  Lewy'  s  work  serves  as  a  starting  point  from  which  we  ulti- 
mately obtain  an  asymptotic  expansion  of  g(z)  at  the  point  cusp.     With  suf- 
ficient care  given  to  the  branch  point,    the  asymptotic  expansion  of  the 
point  cusp  is  found  to  be 

-9  ^0 


(j^nCz-z^))^  (IMz-z^)) 


'[£n{-lln(z  -Zq))] 
+  0  ' 

(Iniz-z^))' 


In  terms  of  x  and  y  this  becomes 

a^{x-x^)^     I         _^  2iln(-lln{x-x^))  (^  +  2lln~-)  Aln{-lln(x-x^))f]j^ 

(£n(x-XQ))^ 


The  expansions  of  g(z)  at  the  cusps  constitute  the  principal  contri- 
bution of  this  dissertation. 


1.     PLASMA  SUSPENDED  IN  AN  AXISYMMETRIC  MAGNETIC  FIELD 

The  physical  situation  of  figure  1  is  that  of  a  fully  ionized  gas  sus- 
pended in  a  perfect  vacuum  by  means  of  a  steady  magnetic  field.     This 
field  is  due  to  a  quadrupole  at  infinity.     The  boundary,    or  interface,    sep- 
arating the  vacuum  from  the  plasma  is  an  axially  symmetric  spindle- 
shaped  figure,    of  the  type  commonly  referred  to  as  "cusped  geometry". 
It  is  the  simplest  such  cusped  geometry,    due  solely  to  the  magnetic  vector 
B  and  the  plasma  itself.     That  is,   there  is  no  electric  field  (E  =  D  =  0) 
and  no  charge  distribution  (p    =  0).     In  the  vacuum  we  may  put  the  magnetic 
permeability  m  =  1  and,    away  from  the  infinite  quadrupole,    B  is  derivable 
from  a  scalar  potential  (^  ■     More  explicitly,    we  have 

B  =  -V  (t> 

The  boundary  consists  of  lines  of  force  of  B  so  that,    if  n  is  the  outward 
normal, 

(1.1)  B.,=-^=0 

on  this  boundary. 

Let  z  denote  distance  along  the  axis  of  symmetry  and  r  denote 
distance  radially  outward  from  this  axis.     Then  at  infinity,    i.  e.  ,    near 
the  quadrupole,    we  find  [18] 

2         2 

(1.  2)  (/)  ~    2z     -  r 

Because  of  Maxwell'  s  equations  we  obtain  [1,  17] 

(1.  3)  7  %    =   -    T'-    B  =  0 

throughout  the  vacuum. 
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The  plasma  is  assumed  to  be  perfectly  conducting  fluid  at  rest.     A 
fluid  continuum  at  rest  can  not  support  internal  pressure  gradients.     There- 
fore the  pressure  is  constant  throughout  the  fluid  and  also  at  the  interface. 
Because  the  plasma  is  a  perfect  conductor  any  current  arising  from  ioniza- 
tion will  be  displaced  instantly  to  the  surface  so  that  there  will  be  only  sur- 
face current. 

To  find  the  surface  force  due  to  the  magnetic  field  we  use  Maxwell'  s 
stress  tensor  [1,  13,  17],    written  in  dyadic  form.      For  E  =  0  this  dyadic  is 
[13] 


rs        1 
T 


477- 


BB  -  ^1{B-B) 


where  I   is  the  identity  dyadic.     The  magnetic  force  per  unit  area  on  the 
interface  is 


n    f  =    ' 


4ir 


(B  ■  B) 


n 


where  we  have  used  n  ■  B  =  0,    which  results  from  the  orthogonality  of 
n  and  B.     This  means  that  there  is  magnetic  pressure  which  may  be  ex- 
pressed as 


■F  •  n  = 


(BB) 


Sir 


In  the  steady  state  the  fluid  pressure  p,    which  is  constant,    must  balance 
the  magnetic  force  (-F-  n)  at  the  interface.     This  gives 


=  -F-  n 


(V4>  ■V<t>) 

Stt 


or 


(1.4) 


(^  <j)  ■\7<^)  "  ^""P  "  constant 
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If  the  shape  of  the  boundary  were  known,    our  problem  would 
be  to  solve  Laplace'  s  equation  for  <ji   in  the  vacuum  and  we  would  need 
only  one  boundary  condition.     However,    our  problem  is  also  to  find  the 
shape  of  the  boundary.     Therefore  what  we  have  is  a  free  boundary 
problem,    which  requires  a  second  boundary  condition.     The  full  naath- 
ematical  statement  of  the  problem  thus  asks  for  a  solution  of 

(1.3)  \/^<j)    =  0 
subject  to  the  boundary  conditions 

(1.1)  V<^    •   n  =  0 

and 

(1.4)  (V</>    •  V<^  )  =  Svrp  =  constant 

This  formulation  also  describes  the  free  boundary  problem  of 
a  vapor-filled  cavity  in  steady  flow  of  water.     There  (/>   is  the  velocity 
potential,   p  is  the  vapor  pressure  which  is  constant,    and  the  boundary 
conditions  arise  from  Bernoulli'  s  theorem  and  the  assumption  that 
water  does  not  enter  the  cavity.    [3] 

Were  we  concerned  here  with  the  hydrodynamic  problem,    we 
would  consider  the  possibility  of  restating  the  problem  in  terms  of  the 
stream  function  0  instead  of  the  velocity  potential  (/>  .     Of  course  our 
plasma  problem  is  now  a  mathematical  one  and  we  are  free  to  intro- 
duce ip  simply  as  a  convenient  function.     Therefore  we  define  the  point 
function  ijj  which  satisfies  equation  (1.  3)  identically  by 

90        ^   ^       dip  d(f)        ^    J_       dip 

9r  r        9z  '  9z  r        9r 

As  we  will  soon  see,   ip  is  constant  on  lines  of  force  of  the  field  of  B. 
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Let  n  and  s  denote  distance  in  a  meridian  plane  along  the  out- 
ward normal  and  along  the  surface,    respectively.     Then  (1.  5)  becomes 

_9^       ^^    _9^  9(/)     _     1        d^ 

^'-  ^'  an    "      r       as  as         r       an 

Because  the  curl  of  a  gradient  is  zero  we  have 

2  2 

1       a    i//  \  d   tp  1  dip 

dr  dz  r 

The  boundary  condition  (1.  1)  is 

^    „  ,    -        I  ,  dip    -      dip   -.     .      -I    dip 

0  -  V  0  •  n  = (  -^   n  -  -^  s)  .    n  =  — 7   ^— 

r        as  an  r      ds 


This  equation  holds  for  any  line  of  force  of  the  field  of  B,    which 
tells  us  that  ip  is  constant  on  each  such  line.     In  particular 

(1.  8)  i//  =  constant 

on  the  boundary.     Substituting  (1.  6)  into  (1.  4)  we  obtain 

constant  =\/  ^  ■  \7*^ 


^    I  .  dip   -      dip      I  i'  .   dip  ^       dip^. 

-  (-^—  n  --r—  s)>       =  — r-  l(^— )     +  (-^)    J 
r       as  an         (  2        as  on 


r       an 


so  that 


(1.  9)  —   -r —    =  constant 

r      an 


on  the  boundary.     Equations  (1.  2)  and  (1.  5)  tell  us  that,    at  infinity,   ip  is 

2 
of  the  form  zr    .     We  may  sum  up  this  analysis  by  re-expressing  our 


mathematical  problem  as 
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(1.7)  ^,^_i9t=0 

9r  9z 


with  the  boundary  conditions 

(1.  8)  ip   =  constant 

and 


1  ^  I 

(1.9)  —   -r —     =  constant 

r       9n 


This  problem,    which  has  thus  far  defied  solution,    will  not  be 
solved  here.     What  we  propose  to  do  is  to  find  the  asymptotic  shape  of 
the  boundary  in  the  vicinity  of  the  edge  cusp,    where  z  is  arbitrarily 
small,    and  the  point  cusps,    where  r  is  arbitrarily  small.     By  utilizing 
the  smallness  of  one  of  the  two  variables  we  can  discard  certain  terms 
and  simplify  the  problem.     To  make  such  a  procedure  useful  we  first 
convert  our  differential  equation  to  a  complex  integral  equation.     In 
this  manner  the  free  boundary  problem  described  by  a  partial  differ- 
ential equation  and  boundary  conditions  for  a  potential  function  ip  becomes 
a  free  boundary  problem  described  by  an  integral  equation  for  a  complex 
analytic  function  g(z)  satisfying  a  single  boundary  condition.     The  desired 
asymptotic  descriptions  of  the  cusps  result  from  a  delicate  analysis  of 
the  free  boundary  problem  for  g(z). 
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2.     EQUIVALENCE  OF  THE  REAL  DIFFERENTIAL  EQUATION 
AND  A  COMPLEX  INTEGRAL  EQUATION 

Garabedian'  s  Extension  to  the  four -dimensional  space  of  two  complex  variables 

Rather  than  refer  specifically  to  either  the  plasma  or  hydrodynamic 
problem  described  above,    we  shall  simply  refer  to  the  mathematical  pro- 
blem which  represents  both.     The  formulation  in  terms  of  (//  will  be  used, 
but  in  cylindrical  coordinates  x  and  y,    where  the  x-axis  is  the  axis  of 
symmetry  of  the  cusped  figure  and  y  denotes  distance  radially  outward 
from  this  axis.     Then  our  equation  (1.  7)  becomes 

(2.1)  i/j       +  ip        -  —   i[j      =  0 

XX       yy      y      y 

with 

(2.2)  ^-_o  ,  ^^=1 

y       9n 

on  the  boundary.     The  constants  0  and  1  are  chosen  arbitrarily,    but 
without  loss  of  generality,    because  ip  is  an  unsealed  potential  to  which 

an  arbitrary  constant  can  always  be  added.     Here  n  denotes  the  outward 

2 
normal  to  the  boundary.     At  infinity,   i//  is  of  the  form  xy    . 

If  the  condition  at  infinity  is  dropped  and  the  free  boundary  is 

prescribed,    our  problem  becomes  a  Cauchy  problem.     Even  from  this 

simplified  point  of  view,    however,    the  attainment  of  a  complete  solution 

is  difficult.     What  we  propose  to  do  is  expand  ip  asymptotically  in  the 

neighborhood  of  the  cusps,    thereby  getting  asymptotic  expressions  for 

the  cusp  contours.     The  asymptotic  approach  we  wish  to  use  for  small 

x  at  one  cusp  and  small  y  at  the  other  is  not  easily  applied  to  equation 

(2.  1)  because  we  do  not  have  any  way  of  knowing  the  magnitude  of  the 

derivatives  of  ip  for  small  x  or  y.     Moreover  the  differential  equation 

(2.  1)  is  extremely  difficult  as  it  contains  the  term  l/y  which  becomes 

infinite  at  the  point  cusps.     Garabedian  [10]  simplifies  the  problem 
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considerably  by  analytic  continuation  into  the  four -dimensional  space 
of  two  complex  variables  {z,z-).     Here  the  elliptic  equation  (2.  1)  in 
(x,  y)  becomes  a  hyperbolic  equation  in  {z,  z'-).   Garabedian  [10]  then 
replaces  the  difficult  differential  equation  by  a  much  simpler  integral 
equation  found  by  solving  the  Cauchy  problem  in  terms  of  a  Riemann 
function.     This  expresses  the  completely  unknown  function  i//  in  terms 
of  an  analytic  function  g(z)  whose  order  of  magnitude  is  known. 
To  begin  with  let  us  consider  the  new  variables 

z  =  z  +  iy  ,  z*  =  X  -  iy 

where  x  and  y  are  themselves  complex.     Then,    clearly, 


9z  9        dz-^    _      9 


9x  9z     9x  dz-''       9x  9z  9z* 


9  9        9z  9        9z-         .   i     9  9 

+  -;: — :    ^ =  1 


whence 


9y  9z       9y         9z=''      9y  I    9z  dz-'- 


9  1,9.9.  9  1/9  .9, 


9z  2      '  9x  9y  '  9z-  2        9x  9y 


9^  9^  .       9^ 


+   ^       =     4 


a    2  2  9z9z- 

ox  9y 


These  expressions  are  substituted  into  equations  (2.  1)  and  (2.  2)  to 
yield 


(2.  3)  i//      „.    +  „,  ^     ,„■    (i//     -iP   J  =  0 

ZZ--        2(z-z--)       z       z=- 


with 

(2.4)  ^-0         ,       j-^     -^ 

(z-z--)        9n 

on  the  boundary. 
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Note  that  the  characteristics  of  equation  (2.  3)  are  the  com- 
plex surfaces  z  =  constant  and  z=i=  =  constant.  It  is  well  known  that 
the  solution  to  an  elliptic  equation  is  analytic.  Therefore  i//,  as  the 
solution  of  (2.  1)  is  an  analytic  function  of  x  and  y.  We  have  defined 
X  =  -^{z  +  z='0  and  y  =-^.{z  -  z*)  which  are  analytic  functions  of  z  and 
z''=.  Therefore  ip  (z,  z*)  is  an  analytic  function  of  z  and  z*  and  may 
be  continued  analytically  into  a  neighborhood  of  the  boundary  in  the 
space  of  two  complex  variables  z  and  z*. 

Garabedian  [10]  observes  that,   for  certain  regions,    equation 
(2.  3)admits  a  formal  solution  in  integral  form.     This  type  of  solution 
employs  as  kernel  a  function  known  as  the  Riemann  function  [8,  9,  10]. 
The  Riemann  function  [8,  9,  12]  of  our  equation  is 


/o    i^\                   D/         -    4-    4--\      V(z  -  t*)(t  -  z"')    -_/(z 
(2.  5)  R(z,  Z'~;  t,  t-'')  = 7 — F' — 


t)(z>-  -  t-:')\ 
t*)(z=:=  -t)j 


(z 


which  satisfies  (2.  3)as  a  function  of  (z,  z'-'')  and  the  adjoint  of  (2.  3) 
as  a  function  of  (t,  t"=).     In  addition 


(2.6)  4^=  l_p  9R  -R 


dz         2(z-z-)  '  9t  2(t-t*) 


on  the  characteristic  z*  =  constant, 


^^-  ''  az*  "    2(z-z=:<)^  '  9t=i=     "    2(z-t-)^ 


on  the  characteristic  z  =  constant,    and 
(2.  8)  R(z,z=:=;z,z*)  =  1 

The  function  F(w),    given  by 
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TT/    ^       XT/       1            1      1         ,      V        [1  •  3  •  5  •  •  ■  (2k  -  3)]^w^ 
F(w)  =  F{  -  -:^,    -  — ;  1  ;  w)  =  S 


2  '         2  '      '    "'       "  2k,,  ,,2 

k  =  o  2       (k!  ) 


is  the  hypergeometric  series  with  parameters  ( — tt,-—^,   1)  and 
satisfies  the  ordinary  differential  equation 


w(l  -  w)F"(w)    +    F'  (w)     -    ^^    =  0 


Before  we  can  use  the  Riemann  function  to  obtain  a  solu- 
tion of  (2.  3)  in  the  four-dimensional  space  of  two  complex  variables 
(z,  Z-),   we  must  show  that  both  Stokes'   theorem  in  a  special  form, 
and  Riemann'  s  representation  formula  for  the  Cauchy  problem, 
hold  in  this  space.     As  Stokes'   theorem  is  required  in  the  deri- 
vation of  Riemann'  s  formula,    we  shall  obtain  the  required  form  of 
Stokes'    theorem  first. 

Special  case  of  Stokes'   theorem 

Let  G  be  any  smooth  two-dimensional  surface  whose  boundary 
is  r    and  let  u(z,  z*)  and  v(z,  z*)  be  two  complex  analytic  functions  of 
the  two  complex  variables  z  and  z=;=  in  the  four -dimensional  space  de- 
fined by  these  two  complex  variables.     Then,    as  we  shall  prove,   the 
following  version  of  Stokes'   theorem  holds: 

(2.9)       jj   [-^    u(z,  z*)  +  -^   v(z,  z=:=)]dzdz-  =         [u(z,  z-)dz*  -  v(z,  z-)dz] 


The  four -dimensional  space  may  be  described  by  z,  z,  z='=  and  z*. 
Then,    according  to  the  calculus  of  exterior  differential  forms  (cf.    Buck 
[6]  )  we  are  able  to  express  the  contour  integrals  of  u  and  v  as  surface 
integrals  of  du  and  dv, 
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j     [u(z,  z-)dz-  -  v(z,  z*)dz]  "//(<: 


z-)dz='=  -  v(z,  z*)dz]  -  j  I  (dudz*  -  dvdz) 


1 


-Sf' 


(— —  dz  +  -^zr-  dz  +  -—^    dz=.=  +  -=;r  dz- )  dz* 
9z  9z  9z-'-  9z-- 


G 


-If' 


(-—^  dz  +  -r^  dz  +  -r— ,.  dz-  +  -^==r  dz*  )  dz 
9z  9z  9z--  9z-'- 


Because  of  the  conventions 


2  2 

(dz)       =  (dz*)       =    0    ;     dzdz*     =    -dz*dz 


and  the  Cauchy-Riemann  equations  [5] 


9-        0     ;     4^,  =  0    ;     4^     =  0    ;    ^     =  0 


9z  '     9z==  '      9z  '     9z^ 


this  becomes 


(2.  9)   /  /  (-^  u(z,z*)  +  ^   v(z,z*))  dzdz*  =  /    (u(z,  z*)dz*   -  v(z,  z*)dz) 


r. 


which  is  the  desired  theorem.     The  result  is  seen  to  combine  the 
more  usual  form  of  Stokes'    theorem  with  the  Cauchy  integral  theorem. 
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Figure  2.       Schematic  View  of  Contour  in  Four-Dimensional   Space 
of  two  Complex  Voriables. 
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Proof  of  Riemann'  s  representation  theorem 

Consider  a  portion  of  the  intersection  of  the  cusped  spindle  of 
plasma  with  a  meridian  plane  and  denote  the  arc  of  its  boundary  running 
from  z*  to  z  by  F.     For  convenience  we  will  refer  to  z'-''-  and  z  as  B  and 
A  respectively,    as  indicated  in  figure  2.     In  the  four -dimensional  space 
of  two  independent  complex  variables  (z,  z*)  the  characteristic  plane 
z  =  constant  which  passes  through  z  and  the  characteristic  plane 
z'-^  =  constant  which  passes  through  z='=  will  intersect  in  a  single  point 
P(z,  z'-'').     Connect  P  and  A  by  means  of  any  differentiable  curve  lying 
in  the  plane  z  =  constant  and  connect  B  and  P  by  means  of  any  differ- 
entiable curve  lying  in  the  plane  z'-'~  =  constant.     Then  pass  a  smooth 
surface  through  ABPA.     We  will  refer  to  the  closed  curve  ABPA  as 
F  ,    and  to  the  smooth  surface  which  it  bounds  as  G. 

If  we  multiply  equation  (2.  3)  by  w  and  its  adjoint  equation  by 
ip  and  subtract  one  from  the  other,    and  then  integrate  the  result  over 
G,    with  (t,  t")  as  running  variables  of  integration  instead  of  (z,  z*),   we 
get 


(2.10) 


-Jff^m)  -  ^L^.Jjdt  dt.  =JJj^  (^  1^   -  j^^.) 


d     ,        dd/  1  ,, 

("  -^   -  o,^     ^-,v   '/'W  )\  dt  dt* 


9t=;=  at       2(t-t*) 

Because  L(!//)  -  0  and  L   (w)  =  0  this  integral  is  equal  to  zero. 
Set 

aw  1 

drjj  1 

V    =    -W   -;: +  -tt: ;  vU 

at  2(t-t*) 
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Whenever  u  and  v  as  defined  above  are  analytic  in  G  and  on  F    we  can 
appl^  Stokes'    theorem  to (2.  10).     We  already  know  that  ip{t,  t*)  is  an 
analytic  function  of  (t,  t-)  but  we  have  not  defined  u  .     Let  u  =  R(z,  z*;t,  V^) 
be  the  Riemann  function  given  by  equation (2.  5).     Then  to  is  analytic 
except  for  poles  and  branch  points  that  occur  at  t  =  t-S    t  =  z*,    and  t*  =  z. 

We  know  that  at  most  one  point  of  T  touches  the  real  axis  t  =t . 
Suppose  that  B  is  displaced  along  T  to  slightly  above  the  real  axis. 
Then  no  point  of  T  touches  this  axis.     The  point  P  at  which  t  =  t*  =  z  =  z- 
lies  in  the  surface  t  =  t*  as  does  the  axis  of  symmetry.     As  we  choose 
the  paths  AP  and  PB  arbitrarily,    we  will  choose  them  so  they  do  not 
intersect  the  surface  t  =  t='S    except  at  P.     Thus  no  point  of  F,  or  G 
touches  the  real  axis,  and  at  no  point  except  P  do  we  have  t  =t*.     At 
A  we  have  t  =  z  so  that  for  t=i=  =  t  to  equal  z  on  F  is  impossible  without 
crossing  the  axis  of  symmetry.     Similarly,    at  B  we  have  t    =  z*  so  that 
it  is  impossible  to  have  t  =  z*  on  F  without  crossing  the  axis  of  symmetry. 

Because  t*  =  z*  on  PB  and  t  =  z  on  AP  we  do  not  have  t  =  z='=  or 
i'-'-  =  z  on  these  curves  except  where  t  =  t-'=  which  is  only  at  P.     The 
function  u   =  R  is  analytic  and  single -valued  at  P  and  clearly  has  no 
singularities  at  any  other  point  of  F  or  G.     Therefore  R  =  u  can  be 
considered  analytic  on  F    and  G  and  we  may  apply  Stokes'   theorem 
(2.  9)  to  (2. 10)  which  becomes 

(2.  12)  0  =  I   I  1^  (^  ^   -  jiji-^,  *R)  -^(R^-  ^fR)i  dt  dt* 


'''   -  W^)  *«'  *'  *  (R|r    -  2(tV)*'"*' 
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We  may  then  rewrite  (2.  12)  as 


/f 


<^-  ^^>  .  1    <«(^  -  2(tV)  ^>^^  -  ^<^  -  2(tV)  ^^''^ 


S  ^^^  -  2aV)  ^)^^  ^  uT^^il  -  20^)  ^)^^*  = « 


-.J     R(- 

PB  AP 

where  we  have  made  use  of  the  facts  that  dt  =  0  on  AP  and  dt*  =  0  on 
PB  because  these  are  segments  of  characteristics  on  which  t  =  constant 
and  t*  =  constant  respectively.     Moreover,   by  (2.  7)  we  see  that  the  last 
integral  in  (2. 13)  is  zero,    and  because  0  =  0  on  P  the  last  three  terms 
in  the  first  integral  of  (2.  13)  vanish.     The  remaining  integrals  are  all 
defined  because  of  the  analyticity  of  R  and  ;//. 
Consider  the  contour  integral 


^1= 


J    «<^  ■  2(tV)  *>*  =  «**p  -  «*'b-  X*<-f  -  2UTTJ-))* 


PB  PB 

The  last  integral  on  the  right  side  of  this  equation  vanishes  by  (2.  6). 
Because  R  is  de: 
by  (2.  8)  so  that 


Because  R  is  defined  at  B  and  i//  =  0  on  F  we  find  (Ri//)    =0.     At  P,    R  =  1 

B 


(2.14)  I^  =  0(z,  z*) 
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Because  i//  =  0  on  F  it  is  evident  that 


-if. 


(2.15)  0=^1    R(-|idt  +  ^dt=:=) 


or 


i/'^-f-^i/^l^- 


When  (2.  15)  and  (2. 14)  are  substituted  in  (2. 13)  we  obtain 


(2.16)  i//{z.  Z-)  =  -  ^     I    R(-^dt  -  ^  dt-) 


iS' 


The  curve  F  lies  in  the  real  domain  t='=  =  t  so  that  we  may  speak 
meaningfully  of  such  things  as  ds  =   |  dt  |  ,    the  differential  of  arc  length 
along  r.     Now  we  wish  to  let  B,    the  lower  end  of  T,    approach  the  real 
axis.     Here  i//  and  its  derivatives  are  single-valued  and  R  is  identically 
zero.     Therefore  equation  (2.  16)  remains  valid  as  B  approaches  the 
real  axis.     At  any  point  of  T  let  n  be  the  unit  outward  normal  and  s 
be  the  unit  tangent.     Moreover  let  i     and  i     be  unit  vectors  in  the 

J.  ^ 

X  and  y  directions  respectively.    Then 

dx    r  dy     r 

q    =    1         -    li—    1 

ds      1        ds      2 
dy    r     ,       dx    T 


9i//     r  dijj     7 


From  these  relationships  we  obtain 
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dtp       r-r  -         dijj       dy  dtp      dx 


=y^  ^  .  n  =  41  ^y   +  ^ 


9n       ^  8x       ds             9y     ds 

which,    with 

,  dip      dt»  dip        dt     .  _      .    Jip_  _d>^        J^/^      dx  . 

^  9t       |dtl  "      at*       |dt|       ~  '^  ^  9x      ds           ay       ds  ' 


is  substituted  in  (2.16)  to  give 


2i  J    ^  9n 


(2.17)  i//(z,  z=;=)  =^.  I    R^     |dt 


Our  free  boundary  condition  on  F  is 

dip    ^    (t  -t---) 
9n  2i 

When  this  is  put  into  (2.  17)  it  yields 

(2. 18)  «z,  z*)  =  ijv(z-f)(z.:.-,)  F(;^:[;g:;;)  1*1 


This  is  Riemann'  s  representation  formula  for  our  problem.     The 
result,    without  the  detailed  proof,    is  given  in  [10]. 
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Equivalence  of  the  real  and  complex  equations 

Equation  (2. 18)  defines  0  on  a  two-dimensional  surface  through 
r  in  the  four -dimensional  space  of  z  and  z*.     We  wish  to  continue  the 
right  side  of  (2.  18)  analytically  to  specify  ip{z,  z*)  for  all  values  of  z 
and  z='=  in  a  four -dimensional  neighborhood  of  F.     Reference  [11]  proves 
that  the  free  boundary  F  is  analytic  and  that  it  is  described  by  the 
equation 

(2.19)  z  =  g(z) 

where  g(z)  is  analytic.     Then  x  and  y  are  taken  as  analytic  functions 
of  arc  length  s  so  that 

z  =  x(s)  +  iy(s)  ,  g(z)  =  x(s)  -  iy(s) 


and 


dg       dz  dg    ,  dz    2  ,    ,  ,    ,   ,2  ,    ,  ,    >  >2 


along  F.     More  explicitly 

ds  =\|g'  (z)  dz 


on  F.     This  is  substituted  into  (2.  18),    along  with  z-  =  g(z-),    to  give 

z 


g(z*) 


(2.20)  ^(z,  z=;0  =  4j-         I    V(z-g(t)  )(z*-t)g'(t)    Fl)^~'T:\:^}".i   ]    dt 
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For  all  g(z)  of  practical  interest  any  branch  point  of  the  integrand  of 
(2.  20)  will  lie  outside  of  the  region  in  which  the  path  of  integration 
lies.     Then,    for  each  fixed  z=;=  on  F  both  sides  of  (2.  20)  are  analytic 
functions  of  z  identical  on  F  and  therefore  identical  in  a  neighborhood 
of  F.     For  each  fixed  z  both  sides  of  (2.  20)  are  analytic  functions  of 
z'-''~  which  have  the  same  values  for  z*  on  F  and  therefore  have  the 
same  values  for  all  z='=.     Thus  by  analytic  continuation,  formula  (2.  20) 
holds  identically  for  all  values  of  z  and  z". 
In  the  real  plane  z*  =  z 
z 


dt 


g(z) 


For  z     any  point  of  F  the  above  integral  may  be  written 


pz       pz       pgu; 

g(z)  Zq        g(z^) 


In  I  set  w  =  g(t).  Then  z  =  g(z)  =  g(g(z)  )  on  F  and  by  analytic 
continuation  we  can  invert  to  obtain  t  =  g(w)  and  g'  (w)g'  (t)  =  I. 
Substituting  into  I  ,  we  find  I 
and  that  the  simpler  formula 


Substituting  into  I    ,    we  find  1=1.     This  shows  that  (//(z,  z)  =  real 


.<z,  i,  -  Re  fj,  ffur^mrrrwiu  f  (i^^M^ 


is  valid  for  z     on  F.     Note  that  0  =  0  on  F  because 
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there.     For  a  known  analytic  boundary  F  this  equation  generates  solu- 
tions of  the  Cauchy  problem  (2.  3),    (2.  4). 

The  equation  ;//  =  0  describes  an  entire  line  of  force  of  the  magnetic 
field.     Only  a  portion  of  this  line  corresponds  to  F,    but  our  analytic 
expression  should  hold  on  the  entire  line,    although  g(z)  =  z  only  on  F. 
By  pushing  z  beyond  F  on  the  line  of  force  of  which  F  is  part  we  might 
hope  to  determine  g(z)  asymptotically.     The  problem  is  then  defined 
in  terms  of  an  analytic  function  g(z)  appearing  in  an  integral  solution  of 
the  more  difficult  partial  differential  equation  (2.  3).     The  principal 
ingredient  with  which  this  general  integral  solution  has  been  constructed 
is  the  Riemann  function. 

It  is  of  some  interest  to  note  that  in  the  two-dimensional  cusped 
geometry  problem,    the  Riemann  function  is  identically  1  and  the  boundary 
conditions  are 


*  =  0        ,  * 


an 


on  F.     Moreover  we  are  in  the  plane  where  z*  =  z  whence  z*  =  z  so 
that 


J_  'iiw^'ij 


i//(z,  z)=—  ,1       (l-l)ds=^.|       ds  =  0 


z-'- 


as  required  by  the  boundary  conditions. 

Now,    before  continuing  with  our  analysis  of  the  axisymmetric 
problem  we  must  examine  the  behavior  of  the  hypergeometric  function 
which  appears  in  the  integrand  of  (2.  21). 
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3.     HYPERGEOMETRIC  FUNCTION 

As  is  well  known  the  hypergeometric  equation 
.2 


w(l  -  w) — ^  +ic-{a  +  b  +  l)wl 

dw      L  J 


dF 

dw 


abw  =  0 


is  satisfied  by  the  hypergeometric  function 

^,     ,  ,      ^       r(a+n)r(b+n)r(c)  n 

^^^'  ^'"'^^  ^  ^     r{a)r(b)r{c+n)r{n+i)     ^ 

n  =  o 

where  r(n+l)  =  n!    is  the  gamma  function.     The  hypergeometric  function 
is  absolutely  convergent  for   |w|<l,  and  divergent  for   |w|>l[19].     When 
|w|   =1  the  function  is  absolutely  convergent  if  [7,  14,  16,  19] 

Re(c  -  a  -  b)>0 
Consider  the  identity  [18] 

(3.  1)    c  I  c  -  1  -  (2c  -  a-b  -  l)w     F{a,  b;  c;  w)  +  (c  -  a)  (c  -  b)wF(a,  b;  c+1;  w) 

^  =  c(c  -  1)  (1  -  w)F(a,  b;  c  -  1;  w) 

for 

Re(c  -  a  -  b)>l. 

The  left  side  of  this  equation  is  absolutely  convergent  at  z  =  1  for 
Re(c  -  a  -  b)>  0  but  the  right  side  requires  Re(c  -  a  -  b)>l.    Then  all 
three  hypergeometric  functions  are  absolutely  convergent  at  w  =  1 
and, by  Abel'  s  theorem  [19]  their  limits  as  w— 1  are  their  values  at 

w  =  1. 
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By  letting  w—  1  in  (3.  1)  we  obtain  the  identity  [19] 


F(a,  b;  c;  1)  =       ^^"  ^^  ^"^  '^,^ F(a,  b;  c+1;  1) 

c(c  -  a  -  b) 


which  is  iterated  to  give  [19] 

-I— r        (c- a  +  n)(c  -  b  +  n)  ., .         -r-/      i,       ._        i\ 

F(a,b;c;l)  =  lim.TT      (e  H- n)(c  -  a  -  b  +  n)  ^^"^-    F(a,b;c  +  m;l) 

m— 00  n=0  m— 00 


The  first   limit  is  clearly 

r(c)r(c  -  a  -  b) 


r(c  -  a)r(c  -  b) 


and  the  second  limit  is  shown  [19]  to  be  unity.     Therefore,    if 
Re(c  -  a  -  b  -  1)>  0  we  have 


/o    ox  TT,      u        ^^        r(c)r(c  -  a  -  b) 

(3.  2)  F(a,  b;  c;  1)  = 


r(c  -  a)r(c  -  b) 


and,    in  particular 


1        1  -  4 

(3.3)  F(-^,  -^;   1;  1)  -  — 
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4.     THE  CUSPS 

In  a  previous  section  of  this  paper  we  considered  an  axial 
cross-section  of  the  spindle  as  shown  in  figure  3.     This  cross-section 
has  two  axes  of  symmetry  and  consists  of  four  mirror  image  curves 
obtainable  from  each  other  by  reflection  across  the  axes  of  symmetry. 
The  longitudinal  axis  of  symmetry,    which  intersects  the  point  cusps, 
is  the  x-axis  of  the  complex  z -plane,   and  the  other  axis  of  symmetry 
is  parallel  to  the  y-axis  of  the  z -plane.     The  cross -section  may  be 
displaced  in  the  x-direction  at  will,    but  not  in  the  y-direction.     This 
is  because  the  differential  equation  (2.  1)  governing  the  magnetic 
field  in  its  exterior  is  invariant  with  respect  to  x-translation  but  not 
with  respect  to  y-translation. 

The  curve  F    (T    ),    together  with  those  portions  of  the  axes 

2\  ID 

of  symmetry  lying  above  it  and  to  its  right  (left),    forms  a  line  of 

force  L    (L    )  of  the  magnetic  field  in  which  the  plasma  is  suspended. 

By  scaling  the  curve  F    we  can  place  the  edge  cusp  at  y  =  1  on  the 

y-axis.     This  leaves  F     in  the  first  quadrant,    which  will  prove  quite 

helpful  as  the  starting  point  of  our  analysis  of  the  edge  cusp.     Distinct 

advantages  in  calculation  at  the  point  cusp  are  realized  if  this  cusp 

lies  at  the  origin  and  if  the  curve  describing  one  arc  of  it  lies  in  the 

first  quadrant,    as  is  the  case  when  the  point  cusp  of  F      is  placed  at 

B 

the  origin.     This  will  be  the  starting  point  of  our  point  cusp  analysis. 

For  convenience,    in  our  separate  analyses  of  the  point  and  edge  cusps 

we  will  drop  the  subscripts  A  and  B  from  F     and  F_  because  there 

A  B 

will  be  no  danger  of  confusion.     When  stating  final  results  the  dis- 
tinction will  again  be  noted,    and  the  expansion  for  a  single  case  will 
be  given.     We  begin  our  separate  analyses  of  the  two  cusps  with  the 
simpler  of  the  two,   the  edge  cusp. 
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edge  cusp 


Figure  3.         Axial  Cross-Section  of  the  Plasma  Spindle 
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5.     THE  EDGE  CUSP 


Confluence  of  anal}i:ic  boundary  conditions 

Figure  4  depicts  the  edge  cusp  as  it  will  be  analyzed  here. 
As  shown  in  section  2,    our  problem  is  to  make  the  integral  I     defined 
by  formula  (2.  21)  real  on  the  line  of  force,    L,    at  least  in  the  neighbor- 
hood of  the  cusps,    subject  to  the  condition  that  the  analytic  function 
g(z)  appearing  in  I     satisfy  the  equation  g(z)  =  z  on  F.    There  is  no 
difficulty  in  proving  that  I     is  real  on  F  itself,    but  the  situation  is 
quite  different  on  the  y-axis.     The  limits  of  integration  of  I     are  at 
our  disposal.     We  choose  them  to  be  the  point  of  separation  i  and  a 
variable  point  iy  on  the  imaginary  axis,    so  that  the  path  of  integration 
will  lie  on  the  positive  y-axis.     Then  the  problem  becomes  as  follows 


(5.1)  F 


f(iy-t)(-iy-g(t) 


^)V 


OJ   n(iy-g(t))(-i;-.''V(-iy-t)ay-g(t>)g'(t)    dt  =  real 


on  the  y-axis  for  y  ^  1,    and 

(5.2)  g(z)  =  i 

on  F. 

The  only  function  at  our  disposal  is  the  analytic  function  g(z). 
On  the  segment  of  the  positive  y-axis  from  i  to  iy  it  must  be  an  analytic 
function  satisfying  the  integral  equation  (5.  1)  and  on  F  it  must  satisfy 
g(z)  =  z.     Thus  at  the  edge  cusp  z=  i  we  have  a  confluence  of  analytic 
boundary  conditions.    Our  ultimate  goal  in  this  section  is  to  obtain  the 
asymptotic  shape  of  F  in  the  neighborhood  of  the  edge  cusp.     This  is 
defined  by  the  asymptotic  behavior  of  g(z)  in  that  region. 
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point  cusp 


Figure  4.  Edge  Cusp  at  z-\. 
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We  note  that  if  g(t)  =  pure  imaginary  for  t  on  the  imaginary  axis, 
then(iy-t),    (-iy-g(t)),  (iy-g(t)),    and  (-iy  - 1)  are  all  pure  imaginary.     In 
the  integrand  this  makes  the  argument  of  F  real,    so  that  F  itself  is  real. 
Moreover,   because  g(t)  is  near  -i  for  t  near  i,    (-iy  -  t)(iy  -  g(t))  =  real>  0. 
Of  course,    g'  (t)  is  then  real  because  it  is  the  change  in  a  pure  imaginary 
with  respect  to  a  pure  imaginary.     If  g'  (t)  is  negative,    its  square  root  is 
imaginary,    and  thus  the  whole  integrand  of  I  becomes  imaginary.     Com- 
bined with  dt,    which  is  imaginary,   this  makes  I     =  real.     If  we  have 
g'  (t)>  0,    then  I     =  pure  imaginary,    contrary  to  our  requirement. 

On  the  basis  of  the  above  observations  it  is  reasonable  to  attempt 
to  prove  deductively  that  g(t)  is  pure  imaginary  for  t  on  the  imaginary  axis. 
If  I     is  differentiated  with  respect  to  iy  and  the  resultant  equation  is  solved 


for  V  -g'  (iy),    we  obtain 


^^0 


(5.3) 


^r¥W)    -    -—-^ 


V(2y)(y +  ig(iy)) 


+ 


r     r2F-(w)Ny(-iy-t)(.y-g(t))[^^^-g^^^y-y'^     ] 
U.        L  (iy-g(t))  (-iy-t) 


V  (-iy-t)(iy-g(t)) 


(iy-t)(-iy-g(t)) 


"^         '         (iy-g(t))(-iy-t) 
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The  terms  appearing  here  are  the  same  as  those  appearing  in  I   ,    so  that 
the  bracketed  term  in  the  integrand  is  real  for  g(t)  =  pure  imaginary. 
Moreover,   this  bracketed  term  is  clearly  bounded  by  some  real  constant 
M.     Here  we  make  use  of  the  fact  that  F'  (w)  is  real  for  real  w.     Because 

is  known  to  be  real,    and  no  further  information  concerning  it  will 


dy  

be  needed,    we  regard  it  as  a  known  quantity.     If  we  set  V  -g'  (t)  =  u(t) 


we  may  write  (5.  3)  as 


u(iy)      =     Hpdy)     +      \       H(t,g(t),  iy,  u{t)  ){-idt) 


1 


■'  -J '"' 


g(t)     =      -i     -I      [u(t)]^  dt 


The  integrand  H  satisfies  a  Lipschitz  condition 


H(t,  g(t),  iy,  u^(t)  )      -      H(t,  g(t),  iy,  U2(t)  )|  <M|u^(t)      -      uj,i)\ 


We  propose  to  use  this  condition  in  connection  with  the  method  of  suc- 
cessive approximations,    to  show   that,    for  g(t)  pure  imaginary,    a  unique 
solution  to  (5.  3)  exists.     The  uniqueness  of  the  solution  assures  that  g(t) 
is  pure  imaginary. 

Now  let  us  define  the  recurrence  relation 

iy 

u^^^(iy)      =     H^dy)     +1      H(t,g^(t),  iy,  u^(t)  )(-idt) 

i 

ly 


/ 


/-: 


gj^  +  l(iy)      =      -i      -    I       u„  (t)    dt 

i 
u„  =  1  ,  Sn  ^     pure  imaginary 
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Our  key  observation  is  that  if  u     =  real  and  g     =  pure  imaginary,   then  by 
the  above  arguments  u  =  real  and  g     ,  -,  =  pure  imaginary,   too.     Thus 

if  u  ,  g  converge  to  u  and  g,  then  u  and  g  respectively  must  also  be  real 
and  pure  imaginary.  Making  use  of  the  method  of  successive  approxima- 
tions,   based  on  our  Lipschitz  condition,    we  obtain 

I  u^  ^  ^(iy)  -  u(iy)  I  <      I       I  H(t,  g(t),  iy,  u^(t)  )  -  H{t,  g(t),  iy,  u(t)  )  |  ( -idt) 

<   m|u   (t)  -  u(t)|         (-idt)  =  mIu   (t)  -  u(t)|(y-l) 
n  '  (j  '    n 

i 

By  iterating  this  relationship  we  find 

|u    ^.(t)  -  u(t)|<  M^^^y-D^'^^lu       .(t)-u(t)| 
n+1  —  n-k 

<   M^'-V-l)"^'!   l-u(t)| 

Because  1  -  u(t)  is  bounded  and  (y  -  1)  may  easily  be  chosen  so  that 

I  y  -  1 1  <  (1/  M)  we  conclude  that  u     converges  to  u.     Every  u     =  real 

which  implies  that  every  g     =  pure  imaginary  and  that  the  limit  of 

the  u     is  real.     Therefore  u  =  V    -g'  (iy)  =  real.     This  means  that 
n 

g'  (iy)  =  real  <0,    which  confirms  that  I     =  real  if  and  only  if  g(t) 
is  pure  imaginary  on  the  segment  of  the  imaginary  axis.     The 
argument  of  the  hypergeometric  function  ir 
real  and  bounded  arbitrarily  close  to  zero. 


argument  of  the  hypergeometric  function  in  the  integrand  of  I     is 


Expansion  of  g(z) 

Having  proven  that  g(z)  =  pure  imaginary  on  the  y-axis,  near 
the  edge  cusp,  we  note  the  obvious  fact  that  both  z  +  g(z)  and  z  -  g(z) 
are  pure  imaginary  here.     On  T  the  known  relationship  z  =  g(z)  tells 
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C-plone 


Figure  5.    Mapping  a  Semi-Neighborhood  of  the  Origin  of  the  ^-Plone 
Upon  a  Semi-Neighborhood  of  the  Edge  Cusp. 
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that  z  +  g(z)  =  2x  -  real  and  z-g{z)  =  2iy  =  pure  imaginary  there.     Now 
consider  the  ? -plane  where  ?  =  ?  +  irj.     Let  us  map  a  semineighborhood 
(rj  >  0)  of  the  origin  of  the  ? -plane  conformally  onto  a  semineighborhood 
of  i  to  the  right  of  F  and  the  y-axis  (from  i  to  iy)  so  that  ?  =  0  goes  into 
z  =  i,   the  segment  of  the  C  -  axis  0<  C<  ?'    goes  into  F,    and  the  seg- 
ment of  the  C  -  axis  0>  C  >  ?"  goes  into  the  segment  of  the  y-axis  Kt  <y. 
Then  z  -  g(z)  is  pure  imaginary  on  the  ?  -  axis  and  z  +  g(z)  is  real  on  the 
positive  I  -  axis  and  pure  imaginary  on  the  negative  ?  -  axis. 

Both  -i(z-g(z))  andy~^(z  +  g(z))  are  real  on  the  real  C  -axis. 
Schwarz'  s  reflection  principle  permits  both  to  be  reflected  below 
the  C  -  axis,    whereby  they  become  analytic  in  the  neighborhood  of 
?  =  0.     Consequently  these  functions  may  be  expanded  in  convergent 
power  series 


°°  ,n 


-i(z  -  g(z))      =     Z     A    ?  (A        =      real) 

n  n 

n  =  0 


00 

\/T(z  +  g(z))    =     2    B   ?  (B        =      real) 

n  n 

n=0 


which  may  be  written  more  clearly  as 


(5.4)  z    -    g(z)       =      2     iA    ^"^ 

n^O       " 

(5.5)  z     +    g(z)     =      Z     B   ^'^  "  2 


n=0 


n 


If  we  set 

B  A 

C  .      (-l)"" ^      •      C  =      (-I)''— ^ 

^2n  -  1  ^     '         2  '         2n  ^     '       2 

we  may  combine  (5.  4)  and  (5.  5)  to  obtain  z  and  g(z).     Then  we  have 
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k_ 

z      =      2       (i)  C     5 

k  =  -l  ^ 

k 

°°  k+1      T 

g(z)    =      2       {-if+V   ? 
k  =  -l  ^ 


From  the  geometry  of  the  cusp  we  see  that  z(?)  and  z'  (?) 
are  bounded  in  the  neighborhood  of  z  =  i.     In  order  that  this  remain 

true  as  C  — 0,    we  must  have  C   ,  =  C,,  =0.     The  known  correspondence 

-1+1  ^ 

z(o)  =  i  tells  us  that  C     =  1.     Thus  far  we  have  not  scaled  the  S -plane 
in  relation  to  the  z -plane.     We  do  so  now  by  requiring  C     -  1.     Then 


we  may  write 


2  00         v+i  2 

z  =  i  -  i?       +       C    C        +         2     (i)       C    ? 

^  k=4  ^ 


_3_  _k^ 

2  00  u .  1  2 

g(z)    =  -i  +  i?     +       Cgr      +        2     (-i)       Cj^? 

k=4 


For  the  right  half  of  the  upper  cusp  we  must  have  C  >    0  in  order  that 
r  be  in  the  first  quadrant.     For  the  left  half  of  the  upper  cusp,    we 
require  C     <  0  in  order  that  F  be  in  the  second  quadrant.     For  the 
lower  cusp  we  take  the  conjugate  of  these  expressions,    which  yield 
complete  analytic  descriptions  of  the  edge  cusp.     They  are  conver- 
gent so  that  we  do  not  have  to  be  satisfied  with  asymptotic  expansions. 

3 

2 
The  most  important  term  in  the  above  is  ?     ,    which  determines 

the  nature  of  the  cusp.     The  exponent  is  exactly  the  same  as  is  obtained 

in  the  two-dimensional  case.    [2] 
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We  are,  of  course,  concerned  with  g{z)  as  a  function  of  z  rather 
than  of  ?.  This  expansion  is  obtained  by  eliminating  ?  from  the  formu- 
las for  z  and  g(z),    which  gives 

3  _n_ 

g(z)       =       -z       +     2iC,(z-i)         +     f      a   (z-i) 

n-4 


If  the  cusp  is  located  at  z  =  z.    =   x    +  iy.  instead  of  z  =  i,   this  becomes 

3  n 

~2~  °°  ~2~ 

(5.6)  g{z)       =      -z      +      2iC^(z-z)       +      Z     a   (z-zj 

3  1  .     n  1 

n=4 


Here  again  the  exponent  2  shows  up.     As  the  C     are  undetermined 
there  is  no  advantage  in  listing  formulas  relating  the  a     and  C    • 

The  above  expansion  of  g(z)  substituted  into  z  =  g(z)  provides 
an  exact  analytic  description  of  the  edge  cusp  in  terms  of  the  complex 
variable  z.       It  is  possible  to  expand  this  equation  asymptotically  in 
terms  of  the  real  variables  x  and  y  as  follows 


^         3C  ^ 
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6.       THE  POINT  CUSP:    PART  I 

Singular  nature  of  the  point  cusp 

Figure  6  depicts  the  point  cusp  as  it  will  be  analyzed. 
Section  4  of  this  paper  introducing  the  cusp  analyses  describes 
how  this  figure  is  arrived  at.     Here  T,    together  with  the  negative 
X-axis  and  the  line  parallel  to  the  y-axis,    extending  upward  from 
the  edge  cusp,  constitutes  a  line  of  force  of  the  magnetic  field  in 
which  the  plasma  is  suspended. 

Our  problem  at  the  point  cusp  is  much  like  that  of  the  edge 
cusp.     We  want  the  integral  I_^    defined  by  (2.  21)  to  be  real  on  the 
above  line  of  force,   subject  to  the  condition  that  the  analytic  func- 
tion g(z),    appearing  in  I^.  ,    satisfy  the  equation  g(z)  =  z  on  P. 
Although  I     is  automatically  real  on  P,    a  different  situation  holds 
on  the  X-axis.     As  we  are  free  to  choose  the  limits  of  integration 
of  I   ,    we  take  them  to  be  the  point  of  separation  0  and  a  variable 
point  X  <  0,    so  that  the  path  of  integration  will  lie  on  the  negative 
real  x-axis.     Then  our  requirement  becomes 


r  pr(x-t)(x-g(t)  )  ^ 
J       ^4x-t)(x-g(t)  )^^  ^' 


Iq    =  I    F[;^^  ,;;^_^;/,  ;]V  (x-t)(x-g(t)  )g'  (t)  dt  -real 


Now  clearly 


(x-t)  (x-g(t)  ) 
(x-t)  (x-g(t)  ) 
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Figure  6.      Point  Cusp  at  the  Origin. 
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and  as  was  shown  earlier  in  section  3  on  the  hypergeometric  function 

4 
F(l)      =      -^ 

IT 

because  c-a-b      =      2  >   1.     Thus  we  set 

:    -    ^ 


whereupon  the  mathematical  statement  of  the  problem  becomes  [10] 
(6.1)  I      =     I   \/(x-t)(x-g(t)  )g'(t)    dt      =      real 

on  the  real  negative  x-axis  and 

(6.  2)  g(z)      =      z 

on  r. 

Of  course,    what  we  really  want  is  the  asymptotic  shape  of  T 
in  the  neighborhood  of  the  point  cusp,   which  is  defined  by  the  asymp- 
totic behavior  of  g(z)  in  that  region.     As  in  the  case  of  the  edge  cusp, 
this  anal;>i;ic  function  g(z)  is  the  only  function  at  our  disposal.     It 
must  satisfy  z  =  g(z)  on  T  and  I  =  real  on  the  segment  0>  t  >  x.     At 
the  point  cusp  this  segment  meets  T  so  that  we  have  a  confluence  of 
two  analytic  boundary  conditions. 

The  similarity  of  this  situation  to  that  which  faced  us  at  the 
edge  cusp  might  tempt  us  to  try  the  method  of  solution  that  worked 
there.     However,    differentiating  I  with  respect  to  z    =    x  does  not 
give  us  an  integral  equation  for^  +  g'  (x)  because  our  integral  has 
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square  root  singularities  at  t  =  x  and  at  g(t)  =  x.     Therefore  the  con- 
fluence of  boundary  conditions  at  the  point  cusp  is  in  a  certain  sense 
singular,    and  the  procedure  employed  at  the  edge  cusp  is  not  appli- 
cable. 

Lewy'  s  analysis 

Apropos  the  confluence  of  analytic  boundary  conditions,   the 
best  reference  is  to  Lewy  [15],   who  treats  a  class  of  problems  not 
unrelated  to  ours.     He  considers  the  conformal  mapping  of  the  region 

indicated  in  figure  7,    but  with  the  curve  F  known  and  expressed  as 

2 
y  =  X  f(x),    where  f(x)  is  a  power  series  in  x  convergent  for  sufficiently 

small   |x|  .     This  differs  markedly  from  our  own  problem,    in  which  T 
is  an  unknown  free  boundary  which  we  wish  to  determine.     Lewy  ex- 
amines the  problem  of  mapping  the  upper  half  ? -plane  (?  =  C  +ir]) 
conformally  upon  a  region  of  the  z -plane  so  that  a  semi-neighborhood 
(r7>0)  of  the  origin  of  the  ? -plane  goes  into  a  semi -neighborhood 
(y>  0)  of  the  origin  of  the  z -plane,    g  =  0  goes  into  z  =  0,    a  portion 
0  <  I  <  C    of  the  I -axis  goes  into  T  and  a  portion  ?"<  C  <  0  of  the 
C-axis  goes  into  the  segment  0>  x>  x"  of  the  x-axis.     By  means  of 
difference-differential  equations  and  integral  equations,    Lewy  est- 
ablishes that  this  unknown  mapping  z(?)  can  be  expanded  asymptoti- 
cally,   in  the  neighborhood  of  the  cusp,    in  terms  of  ?  and  (?ln?). 
Thus  he  finds 


°°  m  n 

mn 
m,  n  =  0 


where  m  and  n  are  non-negative  integers  in  his  special  example. 
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Figure  7.      Mapping  a  Semi-Neighborhood  of  the   Origin  of  the 

5-Plane  Upon  a  Semi-Neighborhood  of  the  Point  Cusp. 
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Lewy'  s  theory  ought  to  be  applicable  to  our  problem,    even 
though  we  do  not  know  the  equation  of  T.     In  our  case  the  general 
theory  suggests  that  g(z)  might  have  an  expansion  of  the  form 

00  a  /3 

g(z)      =      S  A        z    ^(z(lnz))    "  (o     >0,/3    >0) 

r,     rnn  m  n 

m,  n,  =0 


This  procedure  was  tried,   but  it  failed,   apparently  because  of  the 
singular  character  of  our  boundary  conditions  at  their  confluence. 
The  essential  additional  complicating  factor  not  present  in  Lewy'  s 
analysis  is  the  presence  in  I  of  a  branch  point.     Even  though  the 
exact  form  of  Lewy'  s  asymptotic  expansion  of  g(z)  does  not  hold 

in  our  case,    it  will  turn  out  that  an  expansion  in  increasing  powers 

]_ 

2 
of  z     ,    ln(-ln  z),    and  (1/ In  z)  is  valid. 

We  shall  start  by  merely  assuming  that  g(z)  can  be  expanded 

in  terms  of  the  form  z    f(z)  for  a     =    real  andf^z)  an  analytic  function. 

For  the  admissible  value  a     =    0  this  would  substitute  one  analytic 

function  f(z)  for  another  one  g(z).     Therefore  no  loss  of  generality 

is  entailed.     Two  specific  facts  which  are  known  about  g(z)  are  that 

g(0)    =    0  and  that  g{z)    =    z  on  r.     The  first  of  these  facts  tells  us 

that  the  leading  term  in  the  expansion  of  g(z)  is  not  a  constant.     The 

relation  g(z)    =    z  on  T  would  not  be  true  if  the  lowest  order  term  in 

g(z)  were  z    {a  i=  I)  or  z(ln  z)     (/3  ^  0).     Therefore  the  leading  term 

of  g(z)  must  be  z,    and  we  may  write 


(6.  3)  g(z)    =    z    +    az'^f(z)  {a  >  1) 
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Branch  point 

Before  proceeding  with  the  evaluation  of  o-   ,    a,    and  f(z)  we 
must  again  turn  our  attention  to  the  integral  equation 

X 


-S' 


I    =1    \i(x-t)(x-g(t)  )g' (t)    dt      =      real 
From  (6.  3)  above  we  know  that 


g' (t)    =    1    +    [aaz^'h(z)    +    az^f'(t)] 


where  the  bracketed  term  is  of  higher  order  than  1,    so  that^  g'  (t) 
may  be  expanded  binomially.     The  remaining  portion  of  the  integrand 
is\J  (t-x)(g(t)-x),    which  clearly  vanishes  when  t  =  x  and  when  g(t)  =  x. 
Because  of  the  square  root  sign  both  x  and  g     (x)  are  branch  points  of 
the  integrand  of  I.     For  convenience  we  give  the  second  of  these  branch 
points  a  name  by  defining 

(6.4)  t=:=     =    g"\x) 

This,    of  course,    satisfies  the  equation  g(t"'=)    =    x.     We  choose  that 
branch  of\/  (t-x)  which  is  real  for  real  t>  x.     As  this  inequality  holds 
throughout  the  path   of  integration  the  branch  point  x  should  not  dis- 
turb us.     We  choose  that  branch  of\|  (g(t)-x)  which  is  real  for  real 
g(t)>  X.     For  real  g(t)<  x  we  swing  around  V^  via  the  upper  half  of 
the  z -plane  so  thatN^  (g(t)-x)  becomes  pure  imaginary.     Clearly  this 
complicates  the  problem  of  making  I  real.     Even  if  x  -  g(t)  is  not 
real  the  effect  will  persist.     Essentially  what  is  involved  is  a  looping 
of  the  path  of  integration  about  the  branch  point.     This  is  analyzed 
in  the  following  heuristic  discussion. 
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n 


2 
Suppose  at  first  that  g(z)  =  z  -a  iz      for  a     =  real  >   0  and 

n    =  integer  >  3.     We  shall  consider  the  variation  of  g(z)  as  z  goes 


from  r  to  the  negative  x-axis.     To  do  so  we  keep  |z|     =  r„  =  constant 
and  swing  z  through  the  upper  half  plane  to  the  negative  x-axis.     For 
ease  of  illustration  we  set  n  =  5.     The  behavior  of  g(z)  for  n  =  5  is 
typical  of  that  for  any  other  integer  value  of  n  >  3.     In  figure  3  we 
have  indicated  five  successive  positions  of  z,    beginning  with  z    on  F 

TT  i  TT  i  TT  i 

A  ■       .u  u  ~^  "3"  "2"^  TTi 

and  nnoving  through  z     =  re    .  z„  =  r.e     ,  z     =  r.e      ,    and  z^  =  r.e 

■iUo(J4U  5U 

5  5     o9i 


2  ifl  2        2 

We  consider  z  and  z      as  vectors  r.e     and  r„    e        ,    indicated  by  the 

0  0  -^ 

solid  lines  in  figure  3.     For  each  z   ,  (k  =  1,  2,  3,  4,  5)  the  long  solid 

line  is  the  path  of  t  from  t  =  0  to  t  =  z,  .     The  dotted  line  in  each  case 

k 

represents  the  path  of  g(t)  from  g(0)  =  0  to  g(z,  )■     Those  points  marked 
t    ,    which  satisfy  g(t,  )  =  z    ,    are  obtained  approximately  from 

K  K  K 

_5_ 


We  note  that  t     is  a  branch  point  of  ^  g(t)  -  z    . 

For  2=2-1  the  paths  of  t  and  g(t)  lie  to  the  right  of  t    as  t  moves 
along  the  solid  line  to  z    and  g(t)  moves  along  the  broken  line  to  g(z  ). 
A  similar  situation  holds  for  z  =  z    •     For  z  =  z     however,    t„  lies  on 

the  path  of  integration  and  the  paths  of  t  and  g(t)  coincide  from  0  to  z    . 
In  order  that  g(t)  reach  g(z„)  precisely  as  t  reaches  z    ,    it  is  neces- 
sary that  g(t)  pass  through  z     before  t  reaches  z    .     For  z  =  z     the 

solid  and  broken  lines  are  again  separate,   but  both  lie  to  the  left  of  t . . 
Thus      g(t)  -  z     has  looped  around  the  branch  point  t   .     This  was 
overlooked  in  the  latter  part  of  reference  [10].     Thus  I  is  not  single 
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Figure  8.     Integrating  Around  Branch  Points. 
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valued  and  does  not  turn  out  to  be  real.     For  z  =  z     =  x<  0,   the  situation 

5 

is  even  worse.     Here  g(x)  doubles  back  from  x.     The  branch  point  t 

o 

lies  on  the  negative  x-axis  below  z^  =  x.     We  note  that  t  is  really  the 
branch  point  of  Vg(t)     -    x  which  we  previously  labelled  t*. 

To  avoid  looping  about  the  branch  point  t     we  must  choose  the 
path  of  t  so  that  it  always  passes  to  the  right  of  t   .     Then  g(t)  will  keep 
to  the  right  of  z    .     This  is  illustrated  in  figure  8.     For  z  =  z     we  cause 

K  o 

the  path  of  t  to  follow  a  half  circle  to  the  right  of  t     so  that  g(t)  will  follow 

a  little  detour  to  the  right  of  z,,.     For  z  =  z,  a  much  greater  deformation 

of  the  path  of  t  is  required  to  keep  to  the  right  of  t    .     This  keeps  g(t)  to 

the  right  of  z  ,. 

Our  interest  is  primarily  in  z_   =  x.     Here  we  lead  the  path  of  in- 

o 

tegration  out  to  and  around  t='=  and  back  uncfer  the  axis  to  x.     This  causes 

g(t)  to  move  out  to  and  around  z     and  back  under  the  axis  to  g(z    ).     We 

would  like  to  be  free  to  use  any  path  of  integration  in  the  upper  half  plane, 

especially  the  real  negative  x-axis  itself.     As  we  have  just  shown  this 

n 

~2~ 
cannot  be  done  for  g(z)  =  z   -  a   iz     .     If  the  branch  point  t*  were  placed 

below  the  real  negative  x-axis,    however,    we  could  choose  the  portion 
of  this  axis  from  0  to  x  as  our  path  of  integration  for  I.     Then  the  integrand 
of  1  would  not  loop  about  t*,  as  real  t  moves  from  the  origin  to  t  =  x.     More- 
over,   there  would  be  no  need  for  leading  the  path  of  integration  back  under 
the  axis,    or  out  around  any  particular  point.     Therefore  we  impose  as  a 
requirement  of  our  analysis  that  the  branch  point  t*  lie  below  the  path  of 
integration  of  I,    which  is  the  negative  real  x-axis. 

In  order  to  implement  this  requirement  we  must  have  some  method 
of  evaluating  t"  directly.     We  know  that  x  =  g(t='')-     It  will  be  convenient  for 
us  to  obtain  the  inverse  t='=  -  h(x),    where  h,    like  g     is  analytic.     We  find 
this  as  follows.     For  z  on  F, 
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(6.  5)  z      =      g(z)  ;  z      =      g(z) 

Here  we  define  g(z)      =      g(z),  so  that  we  have  taken  the  conjugate  of 
the  function  only.     Substituting  one  of  the  relations  (6.5)  into  the  other 
we  obtain 

z      =      g(g(z)  )  ;  z      =      g(g(z)  ) 

These  algebraic  identities  in  z  and  z  tell  us  that  g  is  the  inverse  of  g. 
They  can  be  used  to  find  h(x).     By  definition  x    =     g(V''')    =    g{h(x)  ).     We 
know  that  x   =    g(g(x)  )  because  gg  is  the  identity  function.     Therefore 

h    =    g 

Here,    of  course,   we  substitute  h    =    g  into  t*      =    h{x)  to  obtain 


O!  ~ 

t-      =      g(x)      =      X     +      ax    f(x) 


Now  that  we  have  an  expansion  of  t-'=  we  must  press  on  to  the 
determination  of  more  terms  of  g(z).     This  means  that  we  want  to  find 
a,  a  ,    and  f(z).     Before  a  and  a   can  be  properly  investigated,    some  in- 
formation about  f(z)  must  be  obtained.     Let 


(6.  6)  f(z)      =      f^(z)      +      if2(z) 


where,    for  z     =    x      =      real    >     0,    both  f  (x)  and  f   (x)  are  real,    f  (x)  is 

J.  ^  L 

positive,    and  f  (x)  is  of  lower   order  than    f„(x)  so  that  f(x)  is  real  to 
lowest  order.     These  assumptions  may  imply  that  a  is  multiplied  by  an 
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appropriate  constant.     If  f(z)  contains  any  logarithmic  terms  we  choose 
that  branch  of  each  logarithm  whose  argument  is  0  on  the  positive  real 
X-axis.     We  must  keep  above  the  branch  point.     Therefore  as  z  swings 
over  the  origin  to  the  negative  real  axis,   via  the  upper  half  plane,   the 
imaginary  part  of  In  z  moves  from  zero  to  +  tt  .     We  also  assume  that 
any  power  of  z  which  is  a  factor  of  f(z)  is  absorbed  into  z         Then,    the 
identity  z    =    g(g(z)  )  requires  that  a  be  pure  imaginary,    i.  e.  , 
a    =    +  a^)!.    for  a     =    real  >   0.     Because  T  is  in  the  first  quadrant  and 
z    =    g(z)  for  z  on  F,    it  is  necessary  that  the  minus  sign  be  chosen  for 
a^,    i.e., 


-^0^'    ^0  ^   ° 


We  know  that  a  >   1.     Consider  the  special  case  a     =    \.     Then 


i-'-      =     X     +     a   ixf(x) 


The  distance  from  x  to  t-  is 


It*      -     x|       =     aQ|xf(x) 


where  f(x)  is  not  divisible  by  any  non-zero  power  of  x.     Thus  the  dis- 
tance from  t-  to  X  is  almost  of  the  same  order  as  x,    so  that  the  branch 
point  may  be  extremely  close  to  the  cusp,    a  situation  we  certainly  do  not 
want.     To  exclude  this  possibility  we  will  restrict  a  to 


a>   \ 
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7.     THE  POINT  CUSP:    PART  II 

Subdividing  the  path  of  integration 

Our  procedure  will  be  to  expand  the  integrand  of  I  in  an  asymp- 
totic series,    and  then  to  integrate  the  lowest  order  terms       The  imagi- 
nary part  of  I,    so  obtained,    is  termed  the  error.     By  choosing  exponents 
and  coefficients  in  the  expansion  of  g(z)  in  such  a  way  that  the  biggest 
term  of  the  error  in  I  disappears,    we  determine  the  asymptotic  expan- 
sion of  g(z)  to  lowest  order.     Of  course,    foi-  geometric  reasons,    t*  must 
be  kept  below  the  negative  x-axis  and    we  must  still  have  g(z)     =    z  on  F. 
The  expansion  of  the  integrand  of  I  is  not  easy.     The  term  g'  (t)  is  readily 
expanded  by  the  binomial  theorem,    but  the  expansion  of 


\/(g(t)  -  x)      =      )/t  -  x  +  at"nt) 


a 


is  a  different  matter  entirely.     Because  t-x  is  larger  than  t    f(t)  for 
small  t,    but  smaller  than  t    f(t)  for  t  near  x,    a  binomial  expansion  is 
out  of  the  question.     However,    a  Taylor  series  expansion  about  the 
branch  point  t-'  would  be  acceptable  if  we  retained  enough  terms.     If 
only  a  finite  number  of  terms  are  retained  we  will  find  ourselves  keep- 
ing a  term  in  (t-t='=)  and  discarding  a  term  in  (t-t='=)     for  n  =  finite  >   1, 
although  for  t  close  enough  to  t*  the  (t-t*)  term  becomes  smaller  than 
the  (t-t=i=)     term  is  near  t  =  0. 

The  logical  resolution  of  this  difficulty  is  to  break  the  path  of 
integration  into  two  parts.     If  the  point  of  division  t      is  chosen  appro- 
priately,   it  becomes  possible  to  expandy  g(t)  -  x  binomially  on  the  lower 
section  of  the  path  while  expanding  g(t)  in  a  Taylor  series  on  the  upper 
section  of  the  path.     To  accomplish  the  binomial  expansion  on  the  lower 
path  we  need  |x  -  t     |  >  >  |t       f  (t    )  |  .     To  accomplish  the  Taylor  series 
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Figure  9.      Subdividing  the  Path  of  Integration  at  tp, 
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expansion  on  the  upper  path  we  need  |t=-  -  t     |<<  |x|  .     As  noted  in  section 

6  we  confine  our  attention  to  values  of  a  greater  than  1.     If 

-   ex  ~ 
t    =  Re(t='0    =    X    +    Re(ax    f(x)  )  and  6   is  a  small  positive  quantity  we  may 

I      V^ 
set  t       =  t       +     6  .     This  will  satisfy  our  requirements  if  |  x  |    <<  6   on  the 

II  I     I  ^  ~c 

lower  path  and  6  <  <  |x|   on  the  upper  path.     Therefore  we  set  6     =     |x| 

for  0<e<Q'-l.     As  will  become  evident  later,   it  will  never  be  necessary 

to  specify  e  or  6   further  then  this. 

Now  we  define 

t    +  6 

(7.1)  I^      -      /  \i(t-x)(g(t)  -  x)g'(t)  dt 

0 


(7.2)  I^     =         /      V(t-x)(g(t)  -  x)g'(t)  dt 

t^.6 


Here  I    and  I     are  the  restrictions  of  the  integral  I  to  the  lower  and  upper 
paths  of  integration,    respectively.     Thus  we  may  write 


(7.  3)  I      =      I^     +      1^ 


By  our  choice  of  the  point  t    +  6  we  have  insured  that,   in  the  upper 
interval,    g(t)  can  be  expanded  in  a  Taylor  series  about  the  branch  point 
t='%    as 

g(t)    =    g(t*)    +    g'(t-)(t-t*)    +    ^^^  (t-t^)^    +... 


=    ^      +      2      g         ^-^    (t-t^O" 
n=l         "' 
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Figure  10.       Avoiding  the  Branch  Point 


-  56  - 


From  g(t)  it  is  easy  to  get  g'  (t).     Thus  we  have 

°o     g<")(t*)  n-1 

n=l 


Because  g  is  an  analytic  function,   all  of  its  derivatives  are  defined  at  t*. 
Moreover   |t  -  t-  |   is  very  small  so  that  the  successive  terms  in  (t  -t*) 
shrink  rapidly.     If  the  series  for  g(t)  and  g'  (t)  are  inserted  in  I_,    we 
obtain 


y^^  n=l      ^        '■  n-1 


^ •"  /p\ 

g'(t*)      r\/(t-x)(t-t:=)  jl    +    -j-   ^^-(ill(t-t-)    +    0((t-t*)^) 
t^+6  ^  ^^^'^ 


To  lowest  order  both  g'  (t)  and  the  bracketed  expression  in  I     are  equal  to  1. 
Therefore,    to  lowest  order  I     is 


/^ 


I2     =      I    Mit-xHt-v-)   dt 

t  +6 
1 


-  57 


This  integral  formula  incorporates  the  termV  (t-x)(t-t=''),    which  branches 
at  t  =  t'''.     To  avoid  taking^  (t-x)(t-t-)  through  the  branch  point,    even  when 
t*  =  t  ,   we  construct  a  tiny  semi-circle  (T   )  about  t  ,    lying  entirely  in  the 
upper  half  plane.     For  0  >  t  >  t    the  vector  t-t*  is  as  shown  in  figure  (a) 
above,    and  t-x  is  real  and  non-negative.     Clearly  the  argument  of  t-t* 
lies  between  0  and  tt/  2.     It  is  zero  if  t-   =  t     <  t  and  tt  /  2  if  Re(t)  =  t  . 

For  t    >    t  >  X  the  vector  t-t*  is  as  shown  in  figure  (b)  above,    and 

TT 

t-x  is  real  and  non-negative.     Now  it  is  clear  that-—   <  arg(t-t='=)<  ir  .     Here 
arg(t-t*)  is  7r/2  if  Re(t)  -  t    and  n   if  t-  =  t  >   t.     There,    as',)  (t-x)(t-t*) 
swings  along  the  negative  real  x-axis,    from  0>  t  >  t    along  T     to  t    >t  >  x 


the  argument  of  t-t*  increases  by  an  angle  6  for  which  0  <  0  <  tt  .     Through- 
out this  range  arg(t-x)  -  0.     Then  argi^  (t-x)(t-t='=)  )  increases  by-^  ,    with 

0  <  -y-  £  -^-     The  equal  sign  holds  only  if  t*  =  t  .     In  the  light  of  this  dis- 
cussion we  evaluate  (7.  4)  at  the  limits  of  integration  to  obtain 


2 1     ,2a 


^T-'«^4>-l„,od2,il-'o'""^"""'^''<l''l'"'l«>'>l'> 


■} 


Here,    for  real  x  •-   0,    we  have  made  use  of  the  convenient  shorthand  notation 


f^U)    -    f^(x)      +      if^Cx)      -      f(x) 


where  the  functions  f  (x)  and  f  (x)  are  real  for  x>    0,    but  complex  for  x<  0, 
as  in  (6.  6). 


58  - 


The  lower  portion  of  the  path  of  integration 

Next  we  shall  turn  our  attention  to  I  .     By  our  choice  of  the  point 
t     +     5    we  have  insured  that,    in  the  lower  interval,   the  integrand  of  L 
can  be  expanded  binomially.     First  we  have 


J,,      ,,    ,,,       ,              ,,      J,    ,    at^f(t)          a^t^^f(t)^     ^     Vt^^f(t)^ 
V(t-x)(g(t)-x)        =      (t-x)^l    + y      +     Of 2 

^^^"^^         8(t-x)^  \{i-^r 


^^      .    ^    at"f(t)         a^'"f(t)'     ^         '■'""'-' 

(t-x)    +    — 8(t-x)  "-         ,   .^       . 

(t  -x) 


\  (t-x)    / 


The  other  term  in  the  integrand  is  also  expanded  binomially  to  give 

^i^^)   =  1  +  ^  ^(^"^(*)^  -^'  t^*"^<^)j'  ^  o(t2("-^)f(t)') 


These  two  expansions  are  then  substituted  into  I    to  obtain 


0  0 

0  0 

2  \^^  /       S^"^ 

a^  r   (t-x)t^''"^f(t)f(t) dt  +  0(  J     t^''(t-x)(f'(t))^  dtj 


4 

0  '    0 


Before  this  expansion  of  I.  can  be  evaluated  we  must  impose  the 
conditions  that  tf  (t)  and  tf"(t)  are  of  higher  order  than  f(t)  and  f  (t)  respec- 
tively,   so  that 
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(7    5)  Urn      tf'(t) 


(7.6) 


t-0        f(t) 
lim      tf"(t) 


t-0         f  (t) 


The  first  of  these  relationships  is  a  reasonable  assumption  because  we  have 
already  specified  that  f(t)  contains  no  non-zero  power  of  t.     Therefore  f(t) 
varies  more  slowly  than  t.     Non-existence  of  the  limit 

lim    tf  (t) 


t-0      f(t) 

would  imply  that  f(t)  is  highly  oscillatory,    which  would  be  unnatural  for  the 
smooth  curve  T  we  wish  to  describe  with  g(z)    =    z    +    az    f(z).     From  Lewy'  s 
work  we  suspect  that  f(t)  may  begin  with  sonae  power  of  in  t.     We  will  not 
assume  this,    because  we  wish  to  deduce  it.     The  relationships  (7.  5)  and 
(7.  6)  are,    however,    consistent  with  the  possibility  that  the  leading  term 
of  f(t)  is  logarithmic. 

The  last  assumption  (7.  6),    regarding  f"(t),    can  be  dispensed  with 
because  application  of  (7.  5)  to  1  yields  the  leading  term  of  f(t),    which  is 
enough  to  determine  that  (7.  6)  holds.     However,    we  will  assume  (7.  6)  for  the 
time  being  in  order  to  expand  I  to  several  terms.     Initially,    only  those  terms 
dependent  on  (7.  5)  will  be  utilized.     After  the  first  term  of  f(t)  is  found  and 
(7.  6)  substantiated  we  will  use  the  terms  of  I  dependent  on  (7.  6)  to  obtain 
higher  order  terms  of  f(t). 

With  (7.  5)  at  our  disposal,    it  is  clear  that  the  remainder  term  in  1 
is  smaller  by  order  of  magnitude  than  any  of  the  integrals  preceding  it. 
Moreover,    we  can  estimate  the  magnitude  of  this  term.     Because  of  (7.  6), 
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a  2 

(t    f  (t) )     is  an  increasing  function  of  t  which  takes  its  largest  values  in 

the  neighborhood  of  the  upper  limit  of  integration,    so  that  we  may  write. 


/ 


h*' 


2a 


2a  +2; 


t^'^(t-x)(f(t))^dt  =  o(ixr""if(x)r) 


The  integrands  of  the  first  and  third  integrals  of  I    are  total  derivatives. 
Integration  by  parts  is  used  for  the  fourth  and  fifth  integrals  to  give 


(7.7) 


; 


Y' 


(t-x)t^'^   "^f(t))^  dt 


2a  2a -r.  „ 


v 


1  A^"^'    xt^M 

a  I    2q'+1  2a     j 


f(t)f'(t) 


t  +6 

I 


+ 


a  J  I  2^+1 


h^'    /2an  2a\ 

~-]      ((f'(t))     +f(t)f"(t))    dt 


(7.8) 


t,+6 


2a      ^  /f^^'^^  i^^ 

(t-x)t         'f(t)f'(t)dt    =f^r-T^    -   ^^^;t— 1  f{t)f'(t) 


2a +1 


2a 


-) 


V 


S 


V    /2a+l  2a' 

'       (f'(t)^  +  f(t)f"(t))    dt 


2a +1 


2a 
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Because  of  (7.  5)  and  (7.  6)  the  integrals  on  the  far  right  of  (7.  7)  and  (7.  8) 
are  of  the  order  of  the  remainder  term  in  I  . 

With  the  aid  cf  the  substitution  t  =  x(l  -  t)  the  second  integral  of  I 


1 


becomes 


(7.9) 


f(t)         ,,                 2a   r              20-                     2    dr 
dt     =     -  X        I      (1  -  t)      f(x(l  -  t)  )     


0 


X 


'/ 


u 


1  1 

"rf(x(i-T))'^-x'"j 


0} 


U 


2  N    2c^  m 

f(x(l-T))     S  (m)(-T)     dr 

m=l 


where  N  equals  2a   or  oo  as  2a   is  or  is  not  an  integer,    and 

t  +  6 
w     =    1    -    ( )    =     Ixl'^      "^-  Odxl'^"   |f(x)|).     Then  the  next  to  last 

integral  in  (7.  9)  is  integrated  by  parts  to  obtain 


•/■ 


(7.10)-x2"(   f(x(l-T))'    ^^ 


2n  9 

x      (jfnT)f(x(l-T)) 


to 


1  1 

+  2x^'^^^f  (jenT)f(x(l-T))f'(x(l-T))dT 

U  (J 


x  ""inlxT'-^'^fCx)^  -  2x^'^^\T£nT  -  T+l)f(x(l-T))f' (x(1-t)) 


(J 


+     2x 


2a' +2 


1 

n 


(TjenT-T+l)(f'  (x(1-t))       +   f(x(l-T)f"(x(l-T)))    dx      + 


(J 


■62- 


3a  -e 


+    0(|xr''"^UnIxl|    |f(x)|If'(x)|) 


{a-l-e)x^°'S>n\x\(nx))^  +  2x^'^^^f(x)f' (x) 


"/ 


+      2x^'^"^^  (TfnT-T+l)(f'(x(l-T))^    +   f(x(l-T))f"(x(l-T)))    dT 


u 


3a  -€ 


+    0{|x|  Un|x|   ||f(x)I|f'(x)i) 


The  last  integral  on  the  right  in  (7.  10)  is  of  the  order  of  the  remainder 
term  in  I  .    In  each  integration  we  have  deliberately  chosen  a  constant  of 
integration  which  causes  the  partial  integral  to  vanish  at  the  upper  limit 
of  integration.     This  was  done  to  avoid  such  terms  as  f(0)  and  f  (0)  which 
might  be  infinite,    as  for  instance  if  the  leading  term  of  f(t)  should  turn 
out  to  be  a  logarithm  raised  to  a  positive  power. 

The  second  integral  of  (7.  9)  is  handled  similarly  to  the  first  one. 


(7.11)       -    x^""  I     f(x(l-T))"    (     Z     (    m)(-T)'"l    dT 


I     f(x(l-T))      [    Z     (    mK-T)'") 


2a      N       2  0-  m 

x  2     (  m  )(-l) 

m-1 


^       2 

f(x)  2x 


m+1         m+2 


f{x)f '  (x) 


^  in-,  r  (^w  -  -  -^  ^2 1  (^'  (-(1-^)'  ^  f(xa-))f"(xa-T)))dr  y  + 


u 
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+        0(|x|2^-^|f(x)||f' (x)|)    +    0(|x|'^-'-2^|f(x)|') 


The  last  integral  on  the  right  of  (7.  11)  is  of  the  order  of  the  remainder 
term  in  I  •     Again  we  have  chosen  the  constant  of  integration  to  avoid 
such  terms  as  f(0)  and  f  (0).     Combining  the  results  (7.  7).    (7.  8),  (7.  10), 
and  (7.  11),    we  have 


real  +  -|-  (t-x)t'^f(t) 


t  +6 


0 


2a  ..    .2      I     la-l-e        2a  ..    .2  N  20-    (-1) 
X      f(x)   jfn  X  -X      f(x)     S(m) 


m 


m=l 


m  +  1 


N    2c?  m 

+    2   S  (m)(-l) 

m-1 


/m+S", 
\m+2/ 


f(x)f'(x)y  +    4^  t^''f(t)^ 


t  +6 


^^^^^^-^-^w 


.2a  n 


]■■' 


+    0(|x  f  (x)      )    +    0(   X  f(x)      ) 


Most  of  these  terms  vanish  at  t  =  0.     When  the  above  expression  is 
evaluated  at  the  limits  of  integration  it  becomes 
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2 1     \2a    2air  i 
I^    =    real    -    ^   I  x  | '^ '^f(x)e^^  ^  -  ^^    ^     / |  (c. -1-e)  Un|  x|   |f(x)' 


+      2    (m)i-4TT  f(x)''  +  22  (m)(-l)"'    -^       x     f(x)f' (x) 
,  (m+1)  ,  m+2 

m=l  m=l 


^0  ''''  .,    ,2    2a7r  i      ^qI^I       ""  J  cyf(x)         ,    |x|f'(x)f(x) 

8  '^"^   ^  -  8  ]2(2a-l)  ^^2_^ 


I     i4«-2-2e|  ,2  I     |2a+2|„,.    .|2 

+      0(|x|  f(x)      )    +    0(   X  f  (x)      ) 


Exponents  of  the  lowest  order  terms  in  g(x)  chosen  to  reduce  the  error  in  I 


To  calculate  I  we  need  only  add  I    and  I   .     Before  doing  so,   however, 
we  recall  that  f(x)  =  f,  (x)  +  if„(x)  where  f^(x)  and  f„(x)  are  real  for  x  real 

and  positive,   and  complex  otherwise.     Then  we  may  write 


I     I  20-  -e                                                       2  I     I  2a-  ,    ,2    2Q'7r  i 
I2-    real  +  — ^-^ (fJOO  +  if^) )  e"'^'''  -  — ^ (eUn|x||  +  0(1)) 

I     |2a  -e                                                     2  I     1  20-^,    ,2    2Q'7ri 
ia„  |x|                                                             ^n  1^1      ^^^'   ^ 
l^=    real  -  —^^—2 (f^(x)  +  if2(x) )  e    ^-— ((a-l-e)Un|x 

+0  (1)  ) 


Adding,   we  obtain 
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I     +    I       =    real    +    la    |  x  |  Im(f2(x)e        ) 


a^lxl 


2a    2aiT  i 


r 


(f^(x)^  -  f^{x)^){a-l)   \ln\x\\   +  2if2{x)f^(x)(cv-l-2eX|£n|x 


>w 


+  0(1)1 


The  constant  e  arises  from  the  definition  of  6 


I  a  -e 


X 


This,    in 


turn,    occurs  in  t      +    6  ,    which  is  the  arbitrarily  chosen  point  at  which  the 
path  of  integration  of  I  is  broken  to  define  I    and  I    .     Of  course  I  must  be 

i  ^  r) 

I     I  Za  -e 
independent  of  e  and  5  .     This  means  that  the  terms  including   |x|  and 

{a  -l-2e)  must  vanish  from  I.     As  this  only  happens  for  f   (x)  =  0,    we  put 


f2(x)     =    0(|x|     )  ,  X>0 


Then  I  becomes,    for  a>l, 


I    =    Real 


2 1     \2a    2aTr  i 
"  'f^(x)   (a    -  l)|£n|x| 


+    0(1) 


J 


Because  f, (x)  is  defined  to  be  real  for  real  x  >    0,    and  to  have  no  non-zero 

2 
power  of  X  as  a  factor,    f    (x)  must  be  real  to  lowest  order  for  x    -    real  <  0. 

Then  the  largest  error  term  in  I  vanishes  for  2a    =  integer.     For  a>   I  the 

principal  solution  would  be  o-     =  3/2.     Solutions  for  a    =  2,    5/2,  3,  7/2,    etc. 

are  higher  order  solutions  that  may  be  expected  to  occur  only  in  degenerate 

or  special  examples.     These  higher  values  of  a  correspond  exactly  to  the 

exponents  of  the  higher  order  terms  in  the  expansion  of  g(z)  at  the  edge 
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cusp.     If  C       =    0  at  the  edge  cusp,   these  higher  order  terms  become 
solutions.     We  will  confine  our  attention  to  the  lowest  exponent  a   =3/2. 

Because  f.(x)  is  complex  it  is  convenient  to  our  subsequent  anal- 
ysis to  define,    for  x    =    real  <   0, 


f^^(x)     =    Re  f^(x)       ;        f^^^'^)     =    ^"^  ^i<^) 


Then  I  becomes 


.     2|     ,3 


I    =    Real  +  ^\^/f^^f^2(UnIx|l    -  2£n(f^^) )    -    Mf^f  -  f  J)    +    0(f^^f^2) 


-|-  |x|(f^^fJ2  +  fi2^n)'+  o(|x|^r(x))^) 


To  make  the  largest  error  term  vanish,    we  must  have,  to  lowest  order. 


which  implies 


f,(x)    =    -^ (jfnx  -  27ri) 

1  TT 


For  X  >  0,    f  (x)  is  real  so  that  this  equation  becomes 


-TTfJx)  -7rf-(x)  .^      -1,  2n  . 

12  (£nx     -    27ri)  \/ Un  x)'  +  (2  rr  )2 
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_2 
which  will  hold  if  f-„(x)  contains  the  factor  (In  x    -    tt  i)      .     Therefore, 

for  X  <  0,    we  set,   to  lowest  order 


£n   X 


we  obtain 


(£n  X  -  277  i)  (£n  X  -  27r  i) 


1  (£n  X  -  TT  i)  o        N        2 

£n  X  (Stn  X  -  2  TT  1)  -tt 


1 

2 

TT 

27r^i 
(£n  x)^ 

+ 

(Jfn  x) 

In  X 

(Inx)^ 

Now  we  have  determined  the  leading  term  of  the  asymptotic  expansion 
of  g{z)  -  z  at  the  point  cusp,   namely 


.    3/2 
a^iz 


{In  z) 


This  term  actually  determines  the  shape  of  the  point  cusp.     If,    however, 
we  should  want  higher  order  terms  in  the  expansion,    we  set 
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^11^^^  =  Wnnr  ^  ^i^^^ 


f     (x)    =    5"   +  p    (x) 


where  p^(x)  and  p„(x)  are  real  and  of  higher  order  than  |£n|x||     and 

|jtn|x||      respectively.     Substituting  f(x)    =    f,-,(x)    +    iL_(x)  into  I  we 
obtain  the  following: 


2|     |3    .     ^ 
I     |x|    TT 1     r  2  ^2 

I    =    real    +    5- <   -2  j^n]  jen|x|  1     -    p(x)|jen|x||       +    — 

5Un|x|r   (^ 


a     |x|    -n  i     C  ....  .       .     ..9  Po(x)Un|x|| 

8 


*  °'u 


^*j 


In  general  we  will  find  the  largest  error  term  in  I  to  be 


ITT  a^      X 


„     X  r  .  v_(x)   £n|x|| 

°  -      C       -    v(x)Un|x|r-l    +    -^ 

n  n  1       '       '     " 


8Un|x|r 


TT 


} 


where  n    =    integer,   C     contains  some  power  of  £n  £n  x      and  v, (x)  and  v„(x) 
are  real  functions.     This  error  term  vanishes  for 
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-(C   +D) 

Un   X  J^n  X 


where  D  is  a  constant.     We  determine  D  from 


V  (x)      +      iv   (x)     =    real  function    -    iG(x) 
1  ^ 


Here  iG(x)  is  the  lowest  order  imaginary  term  produced  by  previously 
calculated  terms  of  f  (x).  Then  we  increase  n  by  1  and  repeat  the  pro- 
cedure. 

To  get  all  terms  of  the  asymptotic  expansion,    we  might  attain  a 
certain  elegance  by  assuming  the  "Ansatz" 


-  -1 

g(z)      =      z    +    2     z        S    2       A,        (jen(-jEnz)) 

m  f  mn 

n=3        m=l    {In  z)  11  =  0 


H  m 

where  the  series  in  (j&n{-£n  z) )  /(In  z)     ought  to  be  convergent  and  the 

n 
2 

series  in  z     is  asymptotic.     The  justification  for  the  "Ansatz"  is  that  we 

have  in  fact  carried  out  the  iterative  scheme  described  above,    so  that  we 
know  the  "Ansatz"  holds  for  n  =  3.     Because  of  the  other  possible  values 
of  a  the  remaining  terms  in  the  "Ansatz"  are  plausible.     Whether  or  not 
this  "Ansatz"  holds  rigorously  is  an  unsolved  problem  which  will  not  be 
attacked  here.     In  any  case  we  find 
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9  ^0 

/    X                             ■    '^  I      "  1           2£n(-£n  z)                 4                   4 
g(z)      =      z     -    a  iz    ^  j^ J—    + 

'  ^"^  (£nz)  (fnz)^ 


^    Q,  (fn  (-£n  z)) 


(Jfnz)^ 


real  + 


2  I     |3     .  /■ 

3.  X  I     TT  i  ■ 

-^- -—J-  \    16(£nUn|x||)^    -    v(x)Un|x|| 

sUnUir  I 


V   (x)Un|x||^ 
+    +    0(£nUn|x|l) 


Consideration  of  g(z)  as  determined  above  tells  us  that  we  have  not 

obtained  as  strong  a  result  as  an  asymptotic  series  in  which  the  terms 

3/  2 
contain  successively  higher  powers  of  z.     Nevertheless  we  do  have  z 

multiplied  by  the  convergent  series 


00 


m-1  £ 


2 2       A,      „  (£n(-£n  z)  ) 

m=l     (£nz)  £  =  0 


and  if  our  Ansatz  holds,    the  expansion  of  g(z)  does  contain  successively 

J. 

2 
higher  powers  of  z"  ,    each  multiplied  by  a  convergent  series  in 

H  m 

{ltn{-S.n  z))  /(£n  z)     .     Therefore,   what  we  have  is  an  acceptable  asymptotic 

description  of  g{z).     The  curve  ^j^oi  figure  3  is  given  by  z    =    g(z).     By  re- 
flecting across  the  axes  of  symmetry  of  this  figure  the  other  curves  are 
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obtained.     To  make  these  expressions  compatible  with  the  results  of  the 
edge  cusp  analysis  we  move  the  point  cusp  to  some  arbitrary  point  z    <  0 
on  the  real  negative  axis.     Then  we  have 


-9  ^0 

3/2/  1  2fn(-£n{z-ZQ))  (^+2£n(^)) 

(7.12)      g(z)    =    z-a^Kz-z^)/       l^j-^--^^-  ^  

0  (!ln(z-z^))  (jPn(z-z    )) 


^(£n(-ln{z-z    )))^' 
+  01  " 


(In(z-ZQ))^ 


As  in  the  case  of  the  edge  cusp,   this  complex  equation  can  be  expanded 
directly  in  terms  of  the  real  variables  x  and  y  to  give 


a^{x-x^)^^^   f         _^  2lln(-fln{x-x^))        ^  +  2Sln(-^) 


I    ln(x-: 


""O^  Unix-x^))^  (Inix-x^))^ 


ilni-Hnix-x   )))^' 
+  01 


ilnix-x^))^ 
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8.     SUMMARY  OF  RESULTS 

1.     Proof  of  Riemann'  s  representation  theorem  in  a  four  dimensional  space. 

In  the  four  dimensional  space  of  two  independent  complex  variables 
(z,  z*)  consider  the  equation 


(2.3)  0      ^      +      -TTT-^ rv  i^      -(//,)     =     0 

^zz*  2(z-z*)        z         z* 


with  the  initial  values 


(2.4)  ^      =      0       ;       -j-^    ^      -      1 

(z  -  z--)      9n 


given  on  an  arbitrary  noncharacteristic  curve  V .     The  Riemann  function 


Rfez-t.w  =  \i(^-f')(t-z')  F(<f-;"^.-':') 

(t  -  t--)  (z  -t-)(z--  -  t) 


satisfies  (2.  3)  as  a  function  of  (z,  z*)  and  the  adjoint  of  (2.  3)  as  a  function 
of  (t,t*).     It  also  satisfies 

z  2(z  -  z==) 

on  the  characteristic  plane  z*  =  constant,    and 


R      ..       ~        Tw Tv  R 

z*  2(z-z*) 


on  the  characteristic  plane  z  =  constant,    and 
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R(z,  z='=;  z,  z*)      =      1 

In  section  2  of  this  dissertation  it  has  been  proved  that  the  solution  to 
(2.  3)  and  (2.  4)  above  is  represented  by 


2i   J   ^    an 


i//(z,  Z-)      =      "2-^  I   R-SZ   |dt 

r 


where  n  is  the  outward  normal  to  T  and  |dt|   is  arc  length.     The  result, 
without  the  detailed  proof,    is  given  in  [10]. 

2.     Convergent  series  representation  of  the  edge  cusp 

An  axial  cross  section  of  figure  1  yields  figure  3,    which  is  composed 
of  four  symmetric  curves.     With  the  origin  of  the  z -plane  at  the  center  of 
this  figure  let  any  one  of  these  four  curves  be  F  which  intersects  the 
y-axis  at  z  .     That  portion  of  F  lying  in  an  arbitrarily  small  neighborhood 
of  the  edge  cusp  can  be  written  as  a  convergent  series. 


n 

3/  2      °°  2 

=    g(z)    =     -    z    +    iC(z-z^)^      +2a^(z-z^) 

n=4 


This  was  proved  in  section  5  of  this  dissertation.     The  representation  is 
valid  for  each  of  the  four  curves  constituting  the  upper  and  lower  cusps 
of  figure  3. 

3.     Asymptotic  representation  of  the  point  cusp 

Consider  the  curve  of  paragraph  3  above  which  intersects  the  x-axis 
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at  z    .     In  a  small  neighborhood  of  j.^,   the  curve  T  is  represented 
asymptotically  by 


-            .^                       .            .3/21           -1               ^^"^-^"^--V>          T"^^"^^) 
z     =    g(z)     =    z-a^Kz-z^)/      }j^^^^z^)     -    -— -^    +    -— -2- 

(£n(z  -Zq))  (£n(z  -  z    )  ) 


^Jfn(-jen(z-z    ))) 
+  0" 


2x 


(£n(z-ZQ))^ 


This  was  proved  in  section  7  of  this  dissertation.     The  representation  is 
valid  for  each  half  of  both  the  right  and  left  cusps  of  figure  3. 
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